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Abstract 


A preliminary  study  of  the  hydroelastic  behavior  of  a faired  tow- 
line  consists  of  a literature  survey  pertinent  to  the  dynamics  of  a 
towline  system,  an  examination  of  the  stability  of  a sectional  fairing 
in  the  light  of  the  existing  theories  and  experimental  evidence  of 
hydroelasticity,  and  a simplified  calculation  of  the  viscous  flow  past 
the  AN/SQA-10  Fleet  Type  Towline  with  the  USNUSL  Sectional  Fairing 
(the  mid-section  of  the  fairing  is  DTMB  Hydrofoil  Shape  No.  7)  in  re- 
lation to  the  self-excited  instability  induced  by  the  viscous  wake. 

The  classical  flutter  is  found  to  be  unlikely  for  a section  of  the 
similar  characteristics.  The  self-excited  instability  requires  a more 
detailed  and  rigorous  calculation  of  the  flow  field;  however,  the 
rough  estimation  indicates  that  the  section  can  be  operated  at  a speed 
of  30  knots  normal  to  the  fairing  without  a serious  separation  of  the 
boundary  layer  from  the  fairing.  The  separation  near  the  mid-body  re- 
sulting a thick  viscous  wake  is  believed  to  cause  the  incipient  of 
this  instability.  A further  study  and  a refined  calculation  of  viscous 


flow  field  is  recommended  for  the  prediction  of  separation  as  well  as  a 

initial  step  toward  the  calculation  of  the  drag  coefficient  for  faired 

towlines  at  any  Reynolds  number. 
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dimensionless  distance  between  mid-chord  and  elastic  axis,  in 
half-chord  length 

half-chord  length  in  feet 

chord  length  in  feet 

pressure  coefficient 

thickness  of  body  or  diameter  of  circular  cylinder  in  feet 
frequency  in  cps 

translational  displacement  in  feet 
shape  factor  (=  £*/®) 

amplitude  functions  for  translational  and  rotational  oscillation 
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mass  polar  moment  of  inertia  per  unit  length  about  the  elastic 
axis 

reduced  frequency  (=Wb/U) 

spring  constants  in  translation  and  rotation 

force  and  moment  per  unit  length  at  the  elastic  axis 

mass  per  unit  length  in  slugs/ft 

pressure  in  psf  or  psi 

the  strength  of  source 

Reynolds  number  (=  Ud/y  ) 

Strouhal  number  (=  -jp 
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parameter  or  time  in  seconds 
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velocity  compc  «?nts  due  to  symmetrical  thickness  in  fps 
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weight  per  unit  length  in  lbs/ft 

shape  function  of  symmetrical  thickness  in  feet 

complex  number  (=  x + iy) 

rotational  angle  or  angle  of  attack  in  radians 
shape  factor  (=  $ jg  (3?)*) 

boundary  layer  and  displacement  thicknesses  in  inches  or  feet 

angle  in  degree  or  radians 
momentum  thickness  in  inches  or  feet 

see  Eq.  (3.2) 

constant  or  parameter  (see  Eq.  3.1) 

2 

shear  stress  on  wall  in  lbs/ ft 
density  in  slugs/ft^ 

circular  oscillation  frequency  in  rads/sec 

translational  and  rotational  natural  frequencies  respectively 
= Kh/m,  <*£  = /I*  ) 

mass  density  ratio  (see  Eq.  2.2) 

2 

kinematic  viscosity  in  ft  /sec 
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I.  Introduction 


A faired  towline,  a heavy-duty  cable  enclosed  in  a symmetrically 
streamlined  body,  is  commonly  used  for  towing  naval  underwater  vessels, 
such  as  a Sonar  device  ( 1 J* , and  oceanographic  instrumentations  /~2.  7 
by  a surface  vessel.  This  dynamical  system  consisting  of  the  towing 
vessel,  towline  and  towed  body  is  essentially  subject  to  the  interac- 
tion of  inertial,  hydrodynamic  and  elastic  forces.  The  dynamics  of 
such  a system  is  indeed  a complex  problem  which  in  general  involves  the 
steady  motion,  the  stability  and  the  unsteady  response  under  different 
environment  for  certain  mode  of  operation.  Progress  has  been  made  so 
far  by  one  of  the  following  approximations:  isolating  an  individual 
component  from  the  whole  system,  treating  only  certain  preferred  mode 
of  motion,  and/or  decoupling  the  forces  acting  on  the  system,  Whi cker 
r 3 7 derived  the  equations  of  motion  for  a cable-body  system  in  planar 
motion  in  terms  of  "generalized",  but  empirical,  loading  functions,  A 
mathematical  solution  of  the  unsteady  forced  motion  based  on  the  method 
of  characteristics  was  obtained  to  facilitate  the  discussion  of  the 
propagation  of  disturbance  along  the  cable.  Asymptotic  solutions  for 
the  perturbation  frequencies  were  also  carried  out  on  the  basis  of  the 
small  disturbance  theory.  His  solution  for  the  equilibrium  configura- 
tion in  terms  of  the  loading  functions  under  a steady  towing  is  availa- 
ble for  practical  applications  and  was  recently  extended  by  Clark  /~L  7. 
Strandhagen  / 5 7 derived  the  equations  of  motion  in  a more  general  way 
by  dividing  a towline  into  a finite  number  of  segments  and  no  solution 
was  attempted  for  the  system  of  equations.  In  these  treatments,  only 
the  inertial  and  hydrodynamic  forces  were  considered  and  the  degrees  of 
freedom  of  dynamical  motion  were  also  limited  in  a certain  way.  The 
elastic  force  has  generally  been  neglected  due  to  the  assumption  either 
that  a towcable  is  extremely  flexible  or  that  hydroelastic  stability  is 
a less  serious  problem  for  the  range  of  speed  considered  at  present  as 
compared  with  the  aeroelastic  problem  for  modern  aircrafts.  The  earlier 
work  in  England  for  aircraft-towed  body  and  the  experimental  program  con- 
ducted at  David  Taylor  Model  Basin  have  contributed  significantly  to  the 
understanding  of  some  phase  of  this  dynamical  problem, 

_ 

The  number  in  / / refers  to  the  listed  references, 
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This  preliminary  study  is  focused  on  the  stability  problem  of  a 
towcable  for  high-speed  operation  as  outlined  in  / 6 /•  The  first 
phase  of  this  study  is  a literature  survey  pertinent  to  the  dynamical 
motion  and  the  stability.  The  list  of  bibliography  given  in  Appendix 
I consists  of  four  parts:  (A)  theoretical  and  experimental  work  direct- 
ly concerning  a towcable  system  or  its  components,  (B)  basic  theories 
and  experiments  in  the  related  hydroelastic  problem,  (C)  the  existing 
knowledge  to  predict  or  to  acquire  the  flow  field  around  a faired  or 
unfaired  towline,  towing  vessel,  and  towed  body,  and  (D)  the  mechanics 
of  wakes  which  is  believed  to  be  essential  in  the  study  of  the  stability 
and  the  dynamic  loading  as  well  as  the  unsteady  response  of  the  system. 
The  second  phase  consists  in  an  examination  of  the  flow  field  around  a 
faired  towline  isolated  from  the  system  under  a steady  towing  and  an 
investigation  of  its  hydroelastic  stability  in  the  light  of  the  existing 
theoretical  and  experimental  evidence.  As  far  as  the  flow  field  is  con- 
cerned, the  faired  towline  may  be  treated  essentially  as  a symmetrical 
hydrofoil  with  a very  large  aspect  ratio  and  a variable  sweep  along  the 
cable.  However,  there  are  basic  differences  in  their  elastic  properties. 
Although  the  theory  of  flow  past  a fully-wetted  airfoil  or  hydrofoil  is 
well  developed,  an  analysis  appropriate  to  a faired  towcable  for  the  pur- 
pose of  extending  the  meager  experimental  data  of  the  drag  coefficient 
(including  the  effect  of  the  inclined  angle)  of  a sectional  fairing  to 
a larger  range  of  Reynolds  number  yet  remains  to  be  developed.  The  cav- 
itating  flow  past  a hydrofoil  is  still  in  the  stage  of  active  research. 
The  study  of  the  hydroelastic  stability  inevitably  requires  a detailed 
flow  field  around  the  cable.  As  will  be  discussed  in  details  in  the 
next  section,  the  classical  type  of  hydroelastic  flutter  is  not  likely 
to  occur  in  the  towcable  problem,  except  in  the  presence  of  cavitation; 
however,  the  self-excited  vibration  due  to  the  separation  of  flow  and 
viscous  wakes  at  high-speed  operation  has  to  be  investigated.  Therefore, 
a procedure  of  numerical  computation  for  the  flow  field  around  a faired 
section  normal  to  the  flow  is  suggested  here.  The  potential  flow  is  ob- 
tained by  a direct  method  based  on  the  thin  airfoil  theory  which  pro- 
vides the  pressure  distribution  on  the  surface  of  the  body  for  a given 
shape  of  fairing.  Then  the  viscous  boundary  layer  over  the  surface  is 
solved  for  various  Reynolds  number.  The  separation  of  flow  from  the  body 
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may  be  determined  and  this  condition  is  essential  for  the  study  of  self- 
induced  vibration.  On  the  other  hand,  this  scheme  may  be  either  employ- 
ed to  design  a faired  section  or  to  extend  the  data  of  experimental 
measurement  of  drag  coefficient  outside  the  range  of  Reynolds  numbers 
available  for  the  purpose  of  determining  the  steady  loading  function. 

Due  to  the  short  span  of  time  available  for  this  study  only  the  first 
part  of  numerical  calculation  was  carried  out  for  the  AN/S(*A-10  Fleet 
Type  Towline  with  the  USNUSL  Sectional  Fairing  for  its  pressure  distri- 
bution and  the  implication  to  the  self-excited  instability  will  be  dis- 
cussed. 

Finally,  this  preliminary  study  is  concluded  with  the  suggestion 
for  future  research  on  the  stability  problem  of  a faired  towline  system. 


II . Hydroelastic  Stability  of  a Towcable 

During  the  last  few  years,  considerable  work  has  been  done  to  extend 
aeroelasticity  to  high-density  fluid  or  licuid  and  hydroelasticity  be- 
comes a new  branch  of  naval  science.  A comprehensive  account  of  aero- 
elasticity may  be  found  in  the  books  by  Bisplinghoff , et.al.  /-7  7 and 
by  Fung  /_8  f.  The  hydroelastic  research,  both  theoretical  and  experi- 
mental, has  so  far  been  scattered  around  various  technical  journals. 
Henry,  et.al.  and  Abramson,  et.al.' s papers  / 9,  10  7 which  appeared  in 
the  Journal  of  Ship  Research  in  1959  are  believed  to  be  the  first  pub- 
lished work  on  hydroelasticity.  The  hydroelastic  stability  is  in  general 
concerned  with  the  phenomena  of  divergence  and  flutter.  Divergence  is  a 
static  instability  due  to  the  interaction  of  hydrodynamic  and  elastic 
forces  and  a divergence  speed  is  defined  to  be  the  critical  speed  over 
which  the  hydrodynamic  force  exceeds  the  elastic  restoring  force  for  a 
specific  mode  of  deformation.  For  a faired  towline,  it  is  generally  de- 
signed in  such  a way  that  the  fairing  is  more  or  less  free  to  rotate 
about  the  "bare"  cable  transversely  to  the  free  stream  and  the  sectional 
fairings  provide  a very  weak  restraint  longitudinally  along  the  cable. 
Therefore,  a fully-wetted,  faired  towline  is  always  statically  stable 
for  steady  speed  and  for  zero  speed  provided  that  certain  conditions  for 
the  centers  of  gravity  and  buoyancy  are  satisfied.  This  implies  that 
there  is  no  divergence  for  a towline.  Flutter  is  defined  as  the  dynamic 
instability  due  to  interaction  of  inertial,  hydrodynamic  and  elastic 
forces  and  the  critical  speed  at  which  this  phenomenon  occurs  is  termed 
as  a flutter  speed.  For  a system  linear  in  its  response  to  loading, 
its  stability  to  infinitismal  disturbance  should  provide  the  complete 
characteristics  of  its  flutter  and  the  origin  of  the  force  producing 
this  motion  is  irrelevant.  On  the  other  hand,  for  a non-linear  system 
there  exists  a functional  relation  between  the  forcing  force  and  its 
response.  When  the  magnitude  of  the  force  increases  the  amplitude  of 
motion  it  provides,  the  phenomenon  is  called  self-excited. 

Now  let  us  consider  a symmetrical  two-dimensional  body  moving  at  a 
steady  speed  in  a fluid  at  rest  or  a uniform  stream  past  a symmetrical 
body  which  has  a thickness-to-chord  ratio  between  1 and  0,  i.e.  the 
shape  of  the  body  varies  between  a circular  cylinder  and  a flat  plate 
of  zero  thickness.  The  flow  past  a circular  cylinder  may  be  character- 
ized by  Reynolds  number  as  given  in  the  following  table  / 11  J: 
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Reynold  Number 


Flow  Description 


0-3 
5 - 10 
40 
90 

150  - 300 

300 


Stoke 's  and  Oseen's  flow 

A pair  of  fixed  vortices  of  Flopp's  type 

The  instability  of  vortices 

The  formation  of  Karman's  vortex  street 

The  formation  of  a shear  layer  from  the 

point  of  separation 

The  irregular  and  turbulent  wake 


As  for  the  range  of  Reynold  number  suitable  for  naval  application  a tur- 
bulent wake  is  inavoidable  for  a bluff  body  such  as  a circular  cylinder. 
The  instability  induced  by  the  unsteady  wake  is  generally  described  by 
the  dimensionless  Strouhal  number  defined  as  S = fd/U  where  the  vortex 
shedding  frequency  f may  be  determined  from  S for  a given  Reynold  number. 
The  Strouhal  number  is  also  known  as  the  reduced  frequency  of  oscillation. 
Strouhal  /-12  7 first  showed  in  1878  that  the  Strouhal  number  S = l/6  for 
a cylinder  and  later  Kovasnay  /~ 13_7  and  ftoshko  /~14_7  determined  for 
circular  cylinder  with  viscous  wakes: 


S = 0.212  (1  - — ),  = -7  (2.1) 

where  X = 21.2  for  40  ^ - 150  and  x = 12.7  for  150  F 5000. 

For  higher  Reynolds  number,  the  "body  regime",  the  region  between  the 
point  of  separation  and  fully-developed  turbulent  wake,  responsible  for 
periodic  oscillation  becomes  shorter.  This  "bare"  cable  represents  a 
limiting  case  of  a faired  towline.  The  stability  is  indeed  a serious  pro- 
blem and  the  use  of  bare  cables  should  be  avoided  even  for  a moderate 
speed  of  operation.  The  study  is  outside  the  present  scope.  The  other 
limiting  case  is  a flat  plate  of  zero  thickness,  which  has  a zero  pressure 
gradient  along  the  plate  and  a uniform  velocity  outside  the  boundary  layer. 
An  extensive  account  on  this  problem  has  been  published  in  many  texts, 
e.g,  Schlichting' s text  /~15  7»  A symmetrical  body  of  finite  thickness 
(including  a thin  plate)  with  non-zero  pressure  gradient  along  the  surface 
and  viscous  wake  behind  the  trailing  edge  is  a practical  model  for  present 
study.  The  hydroelastic  instability  or  flutter  of  a thin  symmetrical 
body  in  a steady,  uniform  stream  may  be  classified  as  the  classical  type 
and  the  nonclassical  type.  In  the  analysis  of  the  former  the  hydrodynamic 
force  is  derived  essentially  from  the  clean  potentially  flow,  i.e.  a 
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fully-wetted  surface  without  a flow  separation  and  viscous  wake,  while 
the  latter  involves  the  flow  of  a real  fluid  - the  boundary  layer  and 
viscous  wake. 

(A)  The  Classical  Flutter.  Let  us  examine  the  classical  flutter 
for  a faired  towline.  Following  Henry,  et.al.  /~9  7 a ncn-dijnenrional 
parameter  paramount  in  a flutter  analysis  is  the  sc  called  "mass-density 
ratio",  which  is  defined  as 

4m 

^ = 1Tf>  c2  '2'2> 

where  m is  the  mass  of  the  faired  towline  per  unit  length,  c the  chord 
and  J)  the  density  of  liquid  or  water.  Since  water  is  roughly  one 
thousand  times  denser  than  sea-level  air,  yU.  turns  out  to  be  or  order  of 
1/2  for  hydrofoils  but  50  or  more  for  aircraft  wings.  This  difference 
works  greatly  in  favor  of  submerged  and  fully-wetted  naval  vessels  for 
dynamic  stability.  A summary  of  the  flutter  analysis  is  given  as  follows. 

The  system  may  be  idealized  by  isolating  a section  from  the  towline 
as  a rigid  hydrofoil  supported  by  two  springs  for  the  bending  and  torsional 
displacements,  h(t)  and  Ot (t),  (see  Fig.  1).  The  equations  of  motion  for 
the  system  are  given  by 


mh  + Srt<X  + K^h  = - L, 
+ Io(«  + K «o<  = M. 


(2.3) 


The  flutter  analysis  is  carried  out  by  assuming  the  solutions  for  an  un- 
damped oscillation  (or  the  neutral  stability). 


_ -z  Jl  wt 

OC  — OIq  e 


h = h e' 
o 


1 cut 


(2.4) 


where  the  flutter  frequency  U5  is  a real  number,  and  o<0  and  hQ  are  com- 
plex amplitude  functions.  The  forcing  functions,  L and  M , are  the  hy- 
drodynamic  force  and  moment  for  this  periodic  oscillation.  For  the  clas- 
sical flutter  analysis,  by  assuming  the  small  disturbances  the  forcing 
functions  may  be  expressed  as  a linear  function  of  the  disturbances,  namely 
ot  , h and  their  first  and  second  derivatives,  e.g.  Theordosen  /~16  7 


- 6 - 


* 
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i 


and  many  others  have  shown  that 

L = T b /“h  + U > - bax_7  + 2 rrp  U bC  (k)  /"h  + U x + b(£  - a)‘*7 
and  (2.5) 

M = ^ b2  /“bah  - Ub(£  - a)*  - b2(l/8  + a2)*7 

y 

+ 2TCU  b2  + a 7 C (k)  /“h  + Uo<  + b(^  - a Yj~J , 

H,(2)(k) 

where  k = b/U,  the  reduced  frequency,  C (k)  = — / „-r t-t* — , and 

Hj2;(k)  + iHU)(k) 
l o 

(2)  (2) 

and  IT  are  the  Hankel  functions  of  the  second  kind.  Substituting 

(2.4)  and  (2.5)  into  (2.3),  two  linear  homogeneous  equations  are  ob- 
tained. In  order  to  have  the  undamped  periodic  solutions  assumed  in 

(2.4) ,  the  determinant  formed  by  the  coefficients  of  K and  « has  to  van- 
ish as  follows:  ! 

/•+  1 - -/<(“)2(^)2  yUX  - a - i fl  + (l-2a)C007  - ; 


, / - a + J /“l  + 2a J C(k),  yUi*2  fl  - (^f)2_7  + a2  + J 


- j r\  - a + (2a2  - £)  C(k)7  + (1  + 2a)  ^ 

k - - ^ 


(2.6) 


This  complex  determinant  is  solved  for  the  critical  speed  UF  and  frequency 
at  flutter  with  the  hydrodynamic,  elastic  and  geometric  properties  as 
parameters.  On  the  basis  of  the  dimensional  argument  discussed  before, 
the  mass  density  ratio  u plays  a paramount  role  for  flutter  as  indicated 
here  in  (2.6).  Fig.  2 is  a typical  result  of  the  theory  and  test  data 
from  Woolston  and  Castile  /~17_/.  As  far  as  flutter  is  concerned  the 
sweep  angle  and  a large  aspect  ratio  tend  to  lower  the  minimum  value  of u ; 
however,  the  existing  knowledge  on  hydroelastic  flutter  indicates  that  the 
asymptote  is  located  atyU>0.5.  The  value  of for  the  AN/SQA-10  Fleet 
Type  Towline  with  the  USNUSL  Sectional  Fairing  is  calculated  as  follows. 
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Chord  of  the  fairing,  c = 0.75  ft. 

Salt  water  density,  = 64/g  slugs/ ft^ 

Unit  weight  of  armored  electrical  cable  (SrU*3986A),  w = 2.58  lbs/ft. 

r 

Unit  weight  of  the  fairing  in  air,  w = 5.16  lbs/ft. 

/ b=  4(wr  + Wf)/ pgrrc2  = 4(5.16  + 2.58)/(64)(3.14  U0.75)2  = 0.274. 


For  this  value  of  the  density  ratio,  there  exists  no  possibility  of  the 
classical  flutter  on  the  basis  of  the  present  state  of  the  hydroelastic 
flutter  analysis  and  of  the  experimental  evidence. 


(B)  The  Self-excited  Vibration 

For  a viscous  flow  past  a thin  fairing,  the  existence  of  viscous 
wake  responsible  for  the  self-excited  vibration  /~ 18,  19,  20  7 and  the 
possibility  of  flutter  due  to  cavitation  / 21  7 are  deserved  a further 
study  of  the  flow  field  around  a faired  tovline  at  a wide  range  of 
Reynolds  number  anticipated  in  modern  naval  operation.  In  the  light  of 
the  experimental  studies,  the  implication  to  the  instability  of  a tow- 
line  is  explored  here.  The  part  of  experimental  results  of  /~18,  19,  20  7 
related  to  this  problem  is  summarized  and  is  discussed  as  follows. 


(l)  The  result  of  Ippen,  et.al.  /~18_7: 
The  configuration:  Flate  No, 


— 2"- 


c 


4 


( n 


c 


Yl 


6 


The  thickness-to-chord  ratio,  d/c  = l/l6 

The  approximate  natural  frequency:  45^100  cps  in  water 

90 ~ 200  cps  in  air 

The  maximum  speed  tested:  35  fps  or  20,7  knots. 


Plate 

No. 

Observed 

Max.  Amplitude 
(radians) 

Strouhal 

No.  at 

Max.  Amplitude 

Strouhal 

No.  at 

Max.  Frequency 

Fre  uency 
at  ''lax. 
Amplitude 

(cos) 

4 

0.003 

- 

- 

165 

380 

3 

0.0035 

0.237 

0.251 

S 380 

5 

0.005 

0.212 

0.219 

85,  175,  390 
470,  530,  570,  - 

6 

0.016 

0.146 

0.174 

95,  156,  205 
285,  370,  680 

A schematic  plot  of  the  characteristics  of  self-induced  vibration  is  shown 
in  Fig.  3»  A linear  relation  for  the  Strouhal  number  with  a distribution 
of  discrete  steps  for  the  resonant  response  coindident  with  the  natural 
frequencies  of  the  elastic  body  was  obtained  from  the  experimental  data. 

One  of  the  important  characteristics  related  to  the  instability  problem, 
namely,  the  separation  condition  in  the  neighborhood  of  the  tail  was  not 
carried  out.  However,  the  result  indicated  that  thers  is  no  serious  insta- 
bility within  the  range  of  test  speeds  for  Plate  No,  4,  which  has  a stream- 
lined tail.  The  wake  is  primary  due  to  the  diffusion  of  the  non-separated 
boundary  layers  detaching  from  both  sides  of  the  body  at  the  trailing  edge. 
The  vorticity  in  this  case  is  very  weak;  so  even  if  the  frequency  coincides 
with  the  natural  frequency  of  the  plate,  the  vibration  is  insignificant. 
However,  if  higher  speeds  were  available  in  the  test,  it  is  expected  that 
the  instability  tends  to  have  similar  characteristics  as  Plate  Nos.  3 and 
6 as  the  boundary  layers  are  separated.  Then,  as  the  speed  increases, 
serious  vibration  will  appear  at  the  discrete  natural  frequencies  of  the 
plate  in  the  spectrum  of  the  exciting  frequencies  induced  by  the  vorticity 
in  the  wakes.  Another  interesting  result  should  be  observed  here  that  the 
vibration  for  Plate  No.  5 is  smaller  than  Plate  Nos.  3 and  6 and  possibly 
No.  4 when  the  flow  is  separated.  This  is  due  to  the  fact  that  for  Plate 
No.  5,  the  points  of  separation  are  more  or  less  fixed  while  the  points 
of  separation  for  other  plates  tend  to  oscillate  and  the  oscillations  of 
lift  tends  to  amplify  the  vibration  further.  Should  a fairing  be  designed 
for  a supercavitating  high-speed  flow,  the  trailing  edge  geometry  of  Plate 
No.  5 is  of  particular  interest. 
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(2)  The  result  of  Heskestad,  et.al.  ^~19 
The  test  plate 

UH" 1 


3> 


The  thickness-to-chord  ratio,  d/c  = 1/46. 
The  range  of  test  speed,  10  to  60  fpe. 


K 


Trailing  Edge  Geometry,  0 


180° 

90° 

60° 

45° 

o 

o 

Plate  No. 

1 

6 

7 

8 

9 

Speed  (fps) 

16.5- 

50.4 

12.1- 

59.7 

16.2- 

60.2 

17.1- 

55.3 

54.8 

Observed 

Frequencies 

(cps) 

145- 

580 

100- 
6 00 

160- 

470 

103- 

597 

597 

Strouhal  No. 

0.192 

0.183 

0.158 

0.275 

0.300 

Distance  Between 
Separation  Point 

d 

d 

d 

d 

none 

Amplitude  at  * 
the  Resonance 

100 

190 

380 

43 

none 

Frequency  of  260  cpe 
* Relative  to  100  for  Plate  No.  1. 

A typical  result  is  reproduced  here  in  Fig.  4,  which  is  similar  to  those 
of  Ippen,  et.al.  In  addition,  the  observation  of  the  distance  between 
separation  points  confirmed  that  the  separation  of  the  boundary  layers 
is  essential  prior  to  the  inception  of  this  instability. 

(3)  Gongwer  /~20  J 
The  configuration 

U U 


The  angle  of  sweep,  0^0^30°. 

The  maximum  speed,  45  knots  for  3/8"-plate  and  60  knots  for  3/4"  plate. 
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The  step  behavior  at  resonant  appeared  on  the  curve  for  the  Strouhal 
number  in  the  frequency/velocity  plots  is  again  confirmed  by  the  ex- 
periment. For  a trailing  edge  thickness  less  than  0.07  inch,  no  in- 
stability was  observed  in  the  speed  range.  The  Strouhal  number  for 
the  instability  varies  with  the  thickness  of  the  trailing  edge,  how- 
ever, if  the  width  of  the  viscous  wake  in  the  fully  developed  region 
instead  of  the  trailing  edge  thickness  is  used,  the  Strouhal  number 
is  approximately  constant  with  a value  of  0.18. 

The  result  of  these  experiments  may  be  summarized  as  the  follow- 
ing: 

(1)  The  self-excited  vibration  may  be  characterized  by  the  Strouhal 
number,  which  has  a constant  value  for  a given  body  configura- 
tion and  elastic  properties  as  indicated  by  the  linear  slope  on 
the  frequency  versus  velocity  plot  for  each  plate  under  a given 
restraint  condition. 

(2)  No  definite  limit  at  which  this  instability  is  incipient  has 
been  observed  because  all  experiments  are  more  or  less  empha- 
sized on  the  response  characteristics  of  the  self-excited 
vibration. 

(3)  No  instability  is  observed  for  a thin  plate  with  a streamlined 
tailing  or  a very  small  trailing  edge  angle.  This  evidence 
implies  that  within  the  range  of  speed  tested,  the  flow  is  not 
separated  from  the  body,  then  the  forcing  frequency  by  the  vorti- 
city  in  the  viscous  wake  primarily  due  to  the  diffusion  of  the 
boundary  layers  detaching  from  the  trailing  edge  is  insignificant. 

(4)  The  resonant  vibration  occurs  as  the  natural  frequencies  of  the 
plate  under  a given  restraint  condition  coincide  with  the  fre- 
quency induced  by  the  vorticity  in  the  viscous  wake.  This  behavior 
appears  in  the  step-like  curve  in  frequency/velocity  plots. 

(5)  The  Strouhal  number  may  be  reduced  approximately  to  a constant 
value  for  the  self-excited  vibration  if  the  width  of  viscous  wake 
instead  of  the  thickness  of  body  is  used  in  the  Strouhal  number. 
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i . Flow  Around  Thin,  Faired  Towcables 

Fince  the  classical  type  of  flutter  is  not  likely  to  occur  and 
only  the  incipient  instability  for  the  self-excited  vibration  is 
crucial  for  the  operation  of  a towcable,  an  analysis  of  the  steady 
flow  of  real  fluid  around  a fairing  is  desired  at  this  stage.  .Ir,  ad- 
dition, the  result  of  this  analysis  may  be  used  to  extend  the  drag 
coefficient  presently  available  from  the  test  to  a wider  range  of 
Reynolds  number,  which  is  essential  in  the  study  of  the  e ui librium 
characteristics  of  a steady  towing.  The  steady  flow  has  to  be  solved 
by  an  iterative  scheme  by  incorporating  the  classical  potential  flow 
theory  with  the  boundary  layer  calculation  and  also  isolating  the  ef- 
fect of  sweep  due  to  the  inclination  of  the  towline,  bears  /~22  J 
has  shown  that  for  a straight  cylinder  in  uniform  flow  the  principle 
of  independent  flow  is  valid,  i.e.  the  cross  flow  along  the  cable 
may  be  evaluated  after  the  main  flow  is  solved.  Although  the  inclina- 
tion angle  might  be  much  smaller  than  90°  and  the  curvature  might  net 
be  small,  it  seems  reasonable  to  separate  this  effect  at  this  stage 
of  analysis  and.  treat  it  as  a uniform  flow  past  symmetrical  hydrofoil 
with  a very  large  aspect  ratio. 

The  iterative  analysis  may  be  carried  out  as  follows.  First, 
solve  the  potential  flow  past  a symmetrical  hydrofoil,  from  which  the 
pressure  distribution  along  the  surface  is  derived.  Then,  the  bound- 
ary layer  along  the  surface  is  calculated  for  a given  Reynold  number 
together  with  the  displacement  and  momentum  thickness  and  the  wall 
shear  stress.  Finally,  potential  flow  has  to  be  solved  again  to  ac- 
count for  the  displacement  thickness.  The  final  result  ’will  provide 
the  pressure  distribution,  which  now  dependent  on  Reynold  number,  may 
be  used  to  predict  the  separation  condition  by  recalculating  the 
whole  boundary  layer.  The  drag  coefficient  may  be  obtained  from  the 
shear  stress  distribution  on  the  surface.  Cf  course,  after  this  scheme 
has  been  set  up  the  iteration  may  be  carried  forward  for  further  refine- 
ment. As  a check  of  the  validity  or  accuracy  of  the  analysis,  the 
computed  drag  coefficient  may  be  checked  against  the  drag  coefficient 
obtained  from  experiments  or  field  measurement  in  the  range  of  Reynold 
numbers  available. 
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Potential  Flow  Solution 

For  a thin  symmetrical  hydrofoil,  the  thin  airfoil  theory  or  the 
linearized  potential  theory  (see  Thwaites  /~23 J)  is  suitable  except  the 
troublesome  singular  behavior  at  the  fore  stagnation  point  and  the  trail- 
ing edge.  In  Appendix  II,  the  calculation  of  this  solution  is  carried 
out  for  the  AN/S^-10  Fleet  Type  towline  with  the  USNU5L  Sectional  Fair- 
ing and  the  result  is  given  in  Fig,  5,  However,  a fairing  with  larger 
thickness-to-chord  ratio,  it  is  possible  to  use  a nonlinear  theory  or 
simply  by  a source  distribution  on  the  surface  instead  of  on  the  chord 
line  as  in  the  linearized  theory. 

Laminar  Boundary  Layer 

Thwaites'  method  /~24  7 for  computing  a laminar  boundary  layer  is  a 
simple  and  straightforward  integration  of  the  potential  flow  velocity  on 
the  surface, 

X = = 0.45  U'U~6  f U5dx  (3.1) 


where  x'  is  from  the  fore  stagnation  point  along  the  surface  of  the  fair- 
ing. And  other  characteristics  such  as  St  S*  > H and  xw  may  be 
computed  from  the  functions  provided  by  the  method  in  the  reference.  A 
check  of  laminar  separation  may  be  predicted  by  /” 25  7 as 


■X  min 


i!  dU 
V dx 


(3.2) 


The  relation  for  and  X is  also  given  by  /~24 J, 
The  Transition 


Since  the  critical  Reynolds  number  for  the  instability  of  a laminar 
boundary  layer  has  been  established  by  the  small  perturbation  theory 

Z"25_7»  £ critical  = (_v‘) critical  is  found  to  depend  on  A * There" 

fore,  from  the  vl  found  in  the  laminar  boundary  layer,  the  point  of  in- 
stability, may  be  located  and  then  the  point  of  transition  is  determined 
from  the  given  empirical  relation  based  on  the  rate  of  amplification  of 
turbulent  eddies  (see  /— 15_7) • 

Turbulent  Boundary  Layer 

It  is  recommended  that  Truckenbrodt • s method  [_  25_7  T°r  calculating 
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the  turbulent  boundary  layer  may  be  used.  The  momentum  thickness  0 and 
the  shape  factor  H have  to  be  calculated  simultaneously  starting  from 
the  transition  point  where  0 is  continuous  while  H has  a finite  drop. 

The  separation  may  again  be  detected  by  checking  the  value  of  H = 1.8  to 

2.4. 


IV  Conclusions  and  Recommendations 

On  the  basis  of  the  existing  theory  and  experimental  evidence  in 
hydro elasticity  related  to  the  instability  problem  of  a faired  towline 
under  steady  motion  examined  in  this  preliminary  study,  the  general 
conclusions  are  given  in  the  following: 

(1)  The  classical  flutter  is  not  likely  for  a faired  towline  having  a 
section  similar  to  the  AN/S^A-10  Fleet  Type  Towline  with  the 
USNl'SL  Sectional  Fairing, 

(2)  The  instability  due  to  self-excited  vibration  requires  a further 
investigation  since  the  existing  theoretical  and  experimental  re- 
sults emphasize  the  response  rather  than  the  incipient  problem. 

However,  the  result  provides  the  necessary  conditions  for  the 
self-excited  vibration,  namely,  the  lower  limit  of  Strouhal  number 
(or  the  Strouhal  number  based  on  the  width  of  the  wake)  and  the 
natural  frequency  of  the  fairing  under  the  practical  restraint 
condition, 

(3)  The  sufficient  condition  for  the  incipient  of  the  self-excited 
instability,  observed  in  those  studies  and  examined  here,  is  the 
separation  of  the  boundary  layers  on  the  surface  of  the  fairing  to 
provide  a strong  viscous  wake  and  the  oscillation  of  the  separation 
points  themselves, 

(4)  On  the  basis  of  an  extremely  simplified  calculation  of  the  characteris- 
tics of  the  turbulent  boundary  layer  for  the  AN/SCiA-10  Fleet  Type  tow- 
line  with  the  USNUSL  Sectional  Fairing  at  a speed  of  30  knots  normal 

to  the  fairing,  it  is  found  that  the  turbulent  boundary  layer  is  not 
separated  from  the  surface  up  to  a point  about  ^ inch  from  the  trail- 
ing edge  due  to  the  singular  behavior  in  the  linearized  potential 
flow  solution.  The  pressure  increases  rapidly  thereafter  while  the 
local  velocity  decreases;  Equation  (J)  indicates  that  the  value  of 
p becomes  rapidly  smaller  than  - 0,06,  However,  even  if  the  flow 
separates  in  the  region  the  width  of  the  viscous  wake  is  expected  to 
be  thin  (the  same  order  of  the  distance  between  separation  point)  and 
therefore  it  is  unlikely  that  this  weak  wake  will  cause  the  self- 
excited  type  vibration, 

(5)  For  a three- fold  increase  of  the  Reynolds  number,  as  seen  from  the  cal- 
culation in  Appendix  II  for  the  prediction  of  separation,  the  conclu- 
sion in  (4)  is  still  valid. 


the  self-excited  vibration  and  should  then  be  incorporated  with 
the  study  of  unsteady  motion. 

The  following  recommendations  are  made: 

(1)  An  iterative  analysis  based  on  the  theory  outlined  should  be  carried 
out  and  then  incorporated  with  a numerical  program  suitable  for  a 
high-speed  digital  computer  in  such  a way  that  for  any  nonanalytic 
shape  of  a fairing  the  flow  field  in  terms  of  the  pressure,  boundary 
layer  characteristics  may  be  calculated  for  any  Reynolds  number. 
Then,  this  program  will  serve  as  the  means 

(a)  to  predict  the  possibility  of  separation  at  a desired  Reynolds 
number  and 

(b)  to  calculate  the  drag  coefficient  outside  the  range  of  Reynolds 
number  available  from  the  experiments. 

This  information  is  not  only  essential  to  the  stability  but  also  to 
the  design  or  improvement  of  a fairing  at  high-speed  operation  and 
to  solve  the  equilibrium  configuration  under  steady  towing. 

(2)  An  experimental  program  to  determine  the  elastic  properties  of  a 
number  of  the  existing  towcables,  e.g.  the  natural  frequency  under 
various  restraint  conditions. 

(3)  The  research  in  unsteady  flow  past  a fully-wetted  and/or  a super- 
cavitating  symmetrical  hydrofoil  is  of  primary  interest  in  many  as- 
pects of  the  towcable  problem,  e.g.  the  stability  and  response  pro- 
blem during  unsteady  operation. 
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Appendix  II 

Prediction  of  Flow  Separation  from  the  Pressure 
Distribution  for  the  AN/SQA-10  Fleet  Type  Towline 
with  the  USNUSL  Sectional  Fairing 


The  Pressure  Distribution 


For  a thin  symmetrical  hydrofoil,  the  potential  flow  may  be  obtained 
for  a given  shape  from  the  thin  airfoil  theory  as  outlined  in  /” 23  7 • Th® 
linearized  solution  consists  of  the  symmetrical  solution, 

i 


u - iv  = - 
s S * 


U • C 


w<  O 


Z - S 


* ys  ^ ^ ) " d l * 


(A) 


and  the  chamber  solution. 


u - iv  = U(-  - 1) 
C c z 


c /~y c(>)-*7 


cC 


)2d£  —~ 


IB) 


L*  Jo  z ~ ** 

where  the  thickness  function,  y Ci),  and  the  chamber  function,  y (5  ), 

s c 

are  given  for  0 £ £ c,  z = x+  iy,  c the  chord  and  C is  a real  constant. 
For  a symmetrical  hydrofoil  like  the  fairing  with  zero  angle  of  incidence, 
only  the  symmetrical  part  of  the  solution  is  needed.  It  is  equivalent  to 
a source  distribution  on  the  chord  line,  O-x-c  and  y = 0 of  the 
strength, 

(C) 


q (x)  = 2Uyg  (x) 


proportional  to  the  slope  of  the  surface  and  with  the  condition  of  closure 
that 


r 


q (x)  dx  = 0. 


(D) 


Then,  the  pressure  distribution  may  be  derived  as  a coefficient, 

2u 


c = 
P 


P - Poo 


k ? U2 


U 


(t) 


and  the  velocity  on  the  surface  is  U = 1 - i c » 

m p 


Let  x = | c (1  + cos  0)  and  % = \ c (1  + cos  t),  then  on  the  surface  of 
the  fairing, 

• ir  (2U/c)  y (t) 

a r dt 

cos  0 - cos  t 


(F) 


it  yn  (t)  dt 


cos  0 - cos  t 


(G) 
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Finally,  by  witson's  r- ethod  /~26_7,  the  numerical  computation  may  be  car- 
ried out  from  the  following  series: 


u = - l—l  CSC  N-1cot  (^)ty'  - y’  ) 

* c *■  1SPSM-1  a ■*' 

odd  p 


or 


(o  = (£)  CSC  (5?)  N_1cot  (^)(y^p  - y^_p) 


F c 
1 m 


l$p$N-l 
odd  p 


(H) 

(I) 


in  which  N is  the  number  of  equal  dimensions  taken  for  0$x$c,  m the 
mth  point  on  the  chord  line  and  the  slope  of  the  surface  y (0)  ■ ^ ^ 

= - — sin  0 y (x)  for  0 = The  convenience  of  this  method  lies  on 
the  fact  that  any  numerically  describable  surface  with  continuous  slope 
may  be  solved  and  the  coefficients  of  the  coefficients  of  the  series  for 
N = 20  are  given  in  /-23  7 * Even  a fairing  may  be  described  exactly  by 
a formula,  however,  after  the  displacement  thickness  due  to  the  boundary 
layer  is  added  to  the  surface  for  a revised  computation  of  the  pressure 
distribution  it  is  not  practical  to  approximate  it  to  a new  equation  or 
an  infinite  series. 

The  pressure  distribution  for  the  AN/S'^A-IO  Fleet  Type  Towline 
with  the  USNUSL  Sectional  Fairing  is  calculated  as  given  in  Fig,  5,  The 
pressure  distribution  shows  the  characteristics  of  a semicircular  nose 
and  then  an  approximately  zero-gradient  over  the  major  portion  of  the 
flat  mid-body.  The  adverse  pressure  gradient  near  the  tail  is  the  center 
of  interest  as  far  as  the  separation  and  the  incipient  instability  of 
the  self-excited  vibration  are  concerned.  Unfortunately  the  linearized 
solution  does  not  provide  a correct  value  due  to  its  singular  behavior 
at  the  fore  stagnation  point  and  the  trailing  edge  and  the  only  possible 
approximation  can  be  made  for  the  purpose  of  predicting  a separation  is 
by  the  dashed  lines  indicated  in  Fig,  5.  Some  sample  data  of  computa- 
tions are  given  in  Tables  1 and  2, 


The  Separation 

The  formal,  conventional  method  for  prediction  has  to  follow  the 
method  given  in  Section  III,  First,  the  solution  of  the  laminar  boundary 
equation  is  calculated  in  terms  of  the  boundary  layer,  displacement  and 
momentum  thickness  as  well  as  the  shape  factor,  then  the  point  of  transi- 
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tion  is  determined,  and  finally  the  turbulent  boundary  layer  is  solved 
in  a way  similar  to  the  laminar  boundary  layer.  In  the  process  of  cal- 
culation, the  possibility  of  separation  is  checked  both  for  the  laminar 
and  turbulent  boundary  layers.  Due  to  the  short  span  for  this  prelimi- 
nary study,  it  is  not  possible  to  carry  out  the  solution  formally,  A 
rough  estimation  is  carried  out  here  on  the  basis  of  the  above  procedure, 
but  certain  approximation  has  to  be  assumed  in  order  to  use  some  known 
results.  The  following  analysis  is  based  on  Reference  /~25  J as  outlined 
in  Section  III.  Lets  consider  a speed  of  30  knots  or  50.7  feet  per 
second  normal  to  the  AN/SQA-10  Fleet  type  Towline  with  the  L'SNUSL  Sectional 
Fairing.  For  a semicircular  nose  with  the  Reynolds  number,  ^G0 |d  = 5.2x10^, 
it  is  reasonably  expected  that  the  transition  point  is  located  in  the 
neighborhood  of  a distance  equal  to  one  half  of  the  thickness  from  the  fore 
stagnation  point  and  there  is  no  possibility  of  a laminar  separation  ieven 
for  a circular  cylinder  at  this  Reynolds  number  there).  This  part  of  the 
laminar  boundary  layer  occupies  only  about  8%  of  the  chord,  therefore, 
the  analysis  in  the  following  is  assumed  to  have  a turbulent  boundary 
layer  starting  from  the  fore  stagnation  point  with  a zero  pressure  gradi- 
ent joint  to  the  pressure  distribution  of  the  mid-body.  We  further  assume 
that  9 computed  by  the  flat  plate  formula  may  be  carried  on  to  the  region 
of  adverse  pressure  gradient  and  the  separation  is  determined  by 


where 

9 = 0.036  x(-^~)-1/5  for  5xl05</P  < 107  (K) 

and  IF?c  = Uooc/y  = 3.15x10^.  Let  x = 8.5"  chosen  for  the  calculation  and 
it  is  found  that 

r = - 0.045 

Therefore,  it  is  not  expected  the  boundary  layer  to  separate  from  the  sur- 
face at  a station  | inch  from  the  trailing  edge.  In  this  approximation, 

9 is  underestimated  in  the  region  of  adverse  pressure  gradient  and  the 
choice  of  the  station  is  quite  arbitrary.  If  a separation  occurs  at  a sta- 
tion, the  thickness  of  the  wake  will  have  the  same  value  approximately  as 
the  distance  between  the  separation  points.  However,  the  calculation  can- 
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not  be  carried  too  close  to  the  trailing  edge  due  to  the  singular  behavior 
of  the  linearized  potential  flow  solution. 

Since  j&q.  v,iO  is  good  for  [R  = 10^,  it  is  expected  that  the  con- 
clusion derived  here  is  still  valid  if  there  is  a three-fold  increase  in 
the  chord  length  for  a geometrically  similar  section. 


'he  Pressure  Distribution  on  the  'JSNUDL  Sect  ionnl  Fairing 


Table  1.  The  Slope  of  the  Mid  Section  of  the  AN/SI^A-IO  Fleet  T.vte 
Towline  with  the  USNUSL  Sectional  lairing 


No.  of 

.. 

cos 

Vc 

*<- 

X 

yi 

y - y 

y. 

Point  i 

(degrees) 

(in.) 

i+i  i- 

i i 

■2 

0 

180 

-1.0000 

0.0000 

0.0000 

0.0000 

175.5 

-.9969 

.0016 

.0144 

.1463 

1 

171 

-.9877 

.0066 

.0594 

.2926 

.2677 

1.704 

166.5 

- .9724 

.0138 

.1242 

.4140 

2 

162 

-.9511 

.0245 

.2205 

.531 0 

.2158 

1.374 

157.5 

-.9239 

.0381 

.3429 

.6298 

3 

153 

.8910 

.0545 

.4905 

.7043 

.1152 

0.0318 

148.5 

-.8526 

.0737 

.6633 

.7450 

4 

144 

.8090 

.0955 

.8595 

.7500 

.0050 

0.0318 

139.5 

.7604 

.1198 

1.0782 

.7500 

5 

135 

.7071 

.1465 

1.3185 

.7500 

.0000 

0 

130.5 

.6495 

.1753 

1.5777 

.7500 

6 

126 

.5878 

.2061 

1.8549 

.7500 

.0000 

0 

121.5 

-.5225 

.2388 

2.1492 

.7500 

7 

117 

.4540 

.2730 

2.4570 

.7500 

.0000 

0 

112.5 

.3827 

.3087 

2.7783 

.7500 

8 

108 

.3090 

.3455 

3.1095 

.7491 

-.0004 

-.0025 

103.5 

.2335 

.3833 

3.4497 

.7496 

9 

99 

.1564 

.4218 

3.7962 

.7482 

-.0038 

-.0242 

94.5 

-.0785 

.4608 

4.1472 

.7458 

10 

90 

0 

.5000 

4.5000 

.7390 

-.0136 

-.0866 

85.5 

+.0785 

.5393 

4.8537 

.7322 

11 

81 

.1564 

.5782 

5.2038 

.7189 

-.0307 

-.1827 

76.5 

.2335 

.6168 

5.55.2 

.7035 

12 

72 

.3090 

.6545 

5.8905 

.6801 

-.0468 

-.2980 

67.5 

.3827 

.6914 

6.2226 

.6567 

13 

63 

.4540 

.7270 

6.5430 

.6245 

-.0645 

-.4106 

58.5 

.5225 

.7613 

6.8517 

.5922 

14 

54 

.5878 

.7939 

7.1451 

.5522 

-.0800 

-.5093 

49.5 

.6495 

.8248 

7.4232 

.5122 

15 

45 

.7071 

.8536 

7.6824 

.4689 

-.0867 

-.5520 

40.5 

.7604 

.8802 

7.9218 

.4255 

16 

36 

.8090 

.9045 

8.1405 

.3794 

-.0922 

-.5869 

31.5 

.8526 

.9263 

8.3367 

.3333 

17 

27 

.8910 

.9455 

8.5095 

.2858 

-.0950 

-.6047 

22.5 

.9239 

.9619 

8.6571 

.2383 

18 

18 

.9511 

.9756 

8.7804 

.1851 

-.1065 

-.6779 

13.5 

.9724 

.9862 

8.8758 

.1318 

19 

9 

.9877 

.9939 

8.9451 

.0736 

-.1165 

-.7416 

4.5 

.9969 

.9984 

8.9856 

.0153 

20 

0 

1.0000 

1.0000 

9.0000 

0.0000 

x = the  distance  from  the  fore  stagnation  point  in  inches 
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Table  2.  A Sample  Calculation  of  C at  x = 0.0594  inch  from  the 
Kore  Stagnation  PSint  (See  Squat ion  1) 


y i+p 

y 1-P 

y 1+p-  y 1-p 

^1^  cscl 

c 

1.374 

1.374 

.8729 

.0318 

-1374 

1.4058 

.2928 

.0000 

-.0318 

.0318 

.0038 

-.0083 

.0000 

-.0083 

-.0007 

-.0866 

.0083 

-.0949 

-.0056 

-.2979 

.0866 

-.3845 

-.0164 

-.5169 

.2979 

-.8148 

-.0249 

-.5869 

.5169 

-1.1038 

-.0228 

-.6779 

.5869 

-1.2648 

-.0152 

— 

.6779 

-.6779 

-.0026 

' 1.0813 

3.072 


